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Abstract 

We examined rates of N2 fixation from the surface to 2000 m depth in the Eastern Tropical South Pacific (ETSP) during El 
Nino (2010) and La Nina (2011). Replicated vertical profiles performed under oxygen-free conditions show that N2 fixation 
takes place both in euphotic and aphotic waters, with rates reaching 155 to 509 |imol N m~^ d~^ in 2010 and 24±14 to 
1 1 8±87 |imol N m~^ d~^ in 2011. In the aphotic layers, volumetric N2 fixation rates were relatively low 
(<1.00 nmol N d~^), but when integrated over the whole aphotic layer, they accounted for 87-90% of total rates 
(euphotic+aphotic) for the two cruises. Phylogenetic studies performed in microcosms experiments confirm the presence of 
diazotrophs in the deep waters of the Oxygen Minimum Zone (OMZ), which were comprised of non-cyanobacterial 
diazotrophs affiliated with nifH clusters 1 K (predominantly comprised of cx-proteobacteria), 1G (predominantly comprised of 
y-proteobacteria), and 3 (sulfate reducing genera of the 5-proteobacteria and Clostridium spp.. Vibrio spp.). Organic and 
inorganic nutrient addition bioassays revealed that amino acids significantly stimulated N2 fixation in the core of the OMZ at 
all stations tested and as did simple carbohydrates at stations located nearest the coast of Peru/Chile. The episodic supply of 
these substrates from upper layers are hypothesized to explain the observed variability of N2 fixation in the ETSP. 
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Introduction 

The efficiency of oceanic carbon (C) sequestration depends 
upon many factors, among which is the availability of nutrients to 
support phytoplankton growth in the illuminated surface ocean. In 
particular, large amounts of nitrogen (N) are required, as it is an 
essential component of proteins, nucleic acids and other cellular 
constituents. Dissolved N in the form of nitrate (NO 3 ) or 
ammonium (NH4^) is directly usable for growth, but concentra- 
tions of fixed N are low (<1 |imol L~^) and often growth-limiting 
in most of the open ocean euphotic zone [1]. Dinitrogen (N2) gas 
dissolved in seawater, on the other hand, is very abundant in the 
euphotic zone (ca. 450 |imol L~^) and could constitute a nearly 
inexhaustible N source for the marine biota. However, only 
certain prokaryotic 'N2-fixers' (or diazotrophs) are able to use this 
N source since they can break the triple bond between the two N 
atoms of the N2 molecule, and convert it into a usable form (i.e. 
NH3) for assimilation. 

The focus of much recent marine N2 fixation research has been 
on the NO 3 -poor environments of the surface tropical ocean, 
where it may sustain up to 50% of 'new' primary production [2,3]. 
The filamentous cyanobacterium Trichodesmium spp., which is 



widespread in the tropical ocean and has a macroscopic growth 
form [4], may fix from 60 [5] to 80 Tg of N per year [6] . Until the 
last decade, this organism was the focus of the bulk of research as it 
is conspicuous and easily collected [4]. However, since then, 
studies of the abundance and diversity of the ni/H gene required 
for N2 fixation have elucidated the importance of unicellular pico- 
and nano-planktonic cyanobacteria [7,8], extending the geograph- 
ical extent of diazotrophy beyond tropical waters [9], and 
potentially narrowing the gap between direct measurements and 
geochemically-based global marine N fixation rates [5]. These 
molecular tools have also revealed the presence of putative non- 
cyanobacterial diazotrophs (possessing and potentially expressing 
the nifH gene) in diverse aquatic environments [10], including 
surface seawater, hydro thermal vents and lakes [11] and 
references therein). In marine waters, these diazotrophs seem to 
be almost ubiquitous [12], but few studies e.g. [13,14] have 
focused on these non-cyanobacterial diazotrophs, and our 
knowledge of their distribution in the ocean and their biogeo- 
chemical importance for the marine N budget is still very limited. 

The N budget for the global ocean is poorly constrained, with 
some suggestions that sinks (denitrification and anammox) exceed 
sources (N2 fixation) [15]. The high energy and iron (Fe) 
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requirements [16,17] of the N2 fixation reaction have implied that 
this process occurs mainly in the large oligotrophic areas of the 
ocean that are depleted in fixed N, and where fixing N2 gives an 
ecological advantage. This may be particularly the case in areas 
which receive high Fe-rich Saharan dust such as the North 
Atlantic [18], or which are under the influence of terrigenous and 
submarine Fe sources, such as the North Pacific near Hawaii 
[19,20] or the South West Pacific [20,21,22,23]. However, recent 
studies [24,25] have hypothesized that N2 fixation might also be 
associated with denitrified surface waters over oxygen minimum 
zones (OMZs), which have measureable NOs", but are depleted 
in N relative to phosphorus (P). This hypothesis has been recently 
confirmed in the coastal surface waters of the Peruvian-Chilean 
upwelling [26,27] as well as throughout the eastern tropical South 
Pacific Ocean (ETSP) [28], where depth-integrated rates over the 
upper water column were comparable to those found in 
subtropical gyres. niJH sequences recovered from these areas 
within the upper 200 m of the ocean were mostly non- 
cyanobacterial and clustered with known heterotrophic sequences 
[26] . This led us to explore N2 fixation in the aphotic zone of the 
ETSP. 

Previous studies conducted in surface waters of the ETSP 
indicated that N2 fixation was highly variable in space and time, 
with depth-integrated rates varying from 10- to 30-fold between 
cruises performed at the same locations [26,28]. Although the 
activities of heterotrophic diazotrophs might potentially be 
contributing to this high temporal variability, very few studies 
have examined the regulation of N2 fixation by heterotrophic 
bacteria in marine waters. Organic C availability has been 
hypothesized to control marine heterotrophic N2 fixation [29] as a 
consequence of the high energy requirements of the reaction, but, 
to our knowledge, the effect of organic molecules on heterotrophic 
N2 fixation has never been studied in OMZs. 

In this study, we investigated N2 fixation along a transect across 
the ETSP in 2010 and 2011 through temperature, oxygen and 
nutrient gradients. We quantified N2 fixation throughout the 0 to 
2000 m depth range in order to evaluate its potential biogeo- 
chemical impact on the marine N budget, and we conducted 
aphotic nutrient addition bioassays in the core of the aphotic 
OMZ in order to investigate which nutrients might control N2 
fixation in this environment. We also phylogenetically character- 
ized the diazotrophs community composition in the core of the 
OMZ and how it responded to some of the nutrient amendments. 

Methods 

Our research was carried out during two cruises in the ETSP, 
aboard the R/V Atlantis in February and March 2010, and the 
R/V Melville in March and AprU 2011. Experiments were 
performed along a transect that began in northern ChUe and ran 
west along 20°S, from the nutrient-rich waters at 82°W to the 
more oligotrophic and low-NOs" waters at 100°W, and returned 
along 10°S (Fig. 1). No specific permissions were required for these 
locations/activities as both cruises took place in international 
waters. This study did not involve endangered or protected 
species. The coastal waters of this region of the ETSP are 
characterized by a permanent wind-driven upwelling of cool 
nutrient-replete water (Fig. 1), which supports high primary 
productivity and a persistent subsurface OMZ, where O2 
concentrations are low enough to induce the anaerobic processes 
of the N cycle, such as denitrification and anammox [30,31,32]. 
These 02-deficient waters are carried by Eckman transport 
westward beyond the limit of our transect. The ETSP is subjected 
to the inter-annual climactic variability of the El Nino-Southern 



Oscillation (ENSO), which modulates the strength of the 
upwelling. The 2010 cruise took place during an El Nino event 
(Multivariate ENSO index: 1.52) and the 2011 cruise during a La 
Nina event (Mutivariate ENSO index: -1.49) fNOAA Climate 
Diagnostics Center, Wolter and Timlin (1993, 1998); data from 
http:/ / www.esrl.noaa.gov/ psd/ enso/ mei.table.html] . 

Hydrographic and nutrient measurements 

Hydrographic and nutrient measurements were performed at 6 
stations in 2010 and 7 stations in 201 1 (Fig. 1). Vertical profiles of 
temperature, chlorophyll <2, fluorescence and dissolved oxygen 
were obtained using a Seabird 911 plus CTD equipped with a 
model 43 oxygen sensor and a Wetlabs ECO-AFL/FL chlorophyll 
fluorometer. Oxygen values were calibrated by micro-Winkler 
[33]. Seawater samples were collected at selected depths using a 
rosette equipped with 24 12-L Niskin bottles. Samples for 
inorganic nutrient (NO3 and P04^ concentrations) analyses 
were collected in acid- washed 20-mL plastic bottles. Nutrient 
concentrations were determined using standard colorimetric 
techniques [34] on a Bran Luebbe AA3 autoanalyzer. Detection 
limits for the procedures were 0.02 to 0.52 |imol L ^ for NOx 
[nitrite (N02")+N03"] and 0.005 to 0.083 |imol L"^ for P04^". 

Vertical profiles of N2 fixation 

Rates of N2 fixation were measured using the ^^N2 tracer 
method [35]. Water samples were dispensed into acid-leached 
4.5-L polycarbonate bottles. During the 2010 cruise, this work was 
exploratory and unreplicated (except for nutrient addition 
bioassays, see below) measurements were made at 12 to 14 depths 
between the surface and 2000 m at stations 1 , 9 and 1 1 . During 
the 201 1 cruise, samples were collected at stations 1, 5, 6, 7, 9 and 
1 1 in triplicates at 1 2 depths between the surface and 2000 m, 
with a specific focus on O2 gradients. Depths were chosen in order 
to sample the oxycline, at least 3 depths within the core of the 
OMZ, as well as an additional 3 within the second increasing 
oxygen concentrations below the OMZ. Most of these depths were 
located in the aphotic zone. 

On both cruises, specific care was taken to avoid O2 
contamination and to perform incubations under strict oxygen- 
free conditions as described in [36]. Briefly, before each profile, 
the 36 4.5-L bottles were filled with deionized water, then the 
deionized water was flushed with argon and finally filled with 
seawater via tubing into the bottom of the argon-filled bottles to 
minimize gas exchange. Bottles were then closed with septa and 
spiked with 3 mL ^^N2 (99 atom % EURISO-TOP) via a gas-tight 
syringe. Each bottle was shaken 30 times to fragment the ^^N2 
bubble and facilitate its dissolution. Recent work has suggested 
that with this method, there may be incomplete equilibration of 
the added ^^N2 gas bubble with the seawater sample, resulting in a 
dissolved ^^N2 concentration in the sample that is lower than the 
equilibrium value assumed in the calculation of ^^N2 fixation rates 
[37]. This may lead to a potential underestimate of N2 fixation 
rates [38,39]. Therefore, the values given in the present study 
should be considered as minimum estimates (discussed below). 
Botdes were then incubated either in on-deck incubators at 
irradiances specific from the sampling depth using blue screening 
and cooled with circulating surface seawater (photic samples), or in 
dark rooms at 12°C or 5°C depending of the sampling depth. After 
incubation, the triplicate bottles from each depth were filtered 
onto precombusted (4 h at 450°C) 25-mm GF/F filters. Filters 
were stored at — 20°C until the end of the cruise, then dried for 
24 h at 60°C and stored dry until mass spectrometric analysis. 
During the 201 1 cruise, an extra 4.5-L bottle was collected at each 
depth of the profile, spiked with ^^N2 and immediately filtered in 
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Figure 1. Location of stations superimposed upon seawater temperature at 75-m-depth. (A) the 2010 cruise (RA' Atlantis), and (B) the 
201 1 cruise (R/V Melville). Station numbering are identical to that of 2 companion papers [28,42]). 
doi:1 0.1 371 /journal.pone.0081 265.g001 



order to determine the initial background 5^^N in the particulate 
organic N (PON) for calculations of N2 fixation rates. During the 

2010 cruise, the value of 5^^N in air (0.00366) was used as a 
reference value for these calculations, which may introduce a 
potential bias, except at Station 1 where ^^N atom % of the PON 
at depth was available. 

Nutrient addition bioassays in the core of the OMZ 

Nutrient addition bioassays of N2 fixation were performed at 
one single depth in the core of the OMZ (based on O2-GTD 
profiles) at 3 stations (Stations 5, 7 and 1 1, between 140- and 450- 
m depth) during the 2010 cruise and at 6 stations (Stations 1, 5, 6, 
7, 9 and 11, between 320- and 475-m depth-) during the 2011 
cruise. All experiments were performed in triplicate and under 
strict oxygen-free conditions (using the argon flushing method 
described above) to avoid inhibition of N2 fixation by oxygen. 
Immediately after collection, bottles were capped with septa and 
amended with nutrients via syringes. During the 2010 cruise, at 
each of the 3 stations, triplicate bottles were left as unamended 
controls, and a second set of bottles was amended with glucose 
to obtain a final concentration of 10 |imol L \ During the 

2011 cruise, triplicate bottles were left as unamended controls, 
and a second set of triplicate bottles was amended with a mixture 
of three simple carbohydrate substrates (39% glucose, 29% 
acetate and 32% pyruvate, final total concentration of 1 |imol 



carbohydrate L-1) to test the effect of a source of dissolved organic 
C (DOC) on N2 fixation. A third set was amended with a mixture 
of three amino acids as a source of both DOC and dissolved 
organic N (DON) (20% leucine, 23% glutamic acid and 56% 
alanine) to reach a final concentration of 1 |imol amino acids L ^ . 
The proportion of each carbohydrate and amino-acid has been 
chosen in order to add the same quantity of organic C in the two 
treatments (4 |imol L ^). A fourth set was amended with ATP 
(source of dissolved organic P, DOP) to reach a final concentration 
of 1 nmol L \ and a fifth set was amended with 8 |imol L ^ of 
NOs" to test its potential inhibitory effect on heterotrophic N2 
fixation. Bottles were then incubated in a dark cold room at 12°C 
for 24 h in order to leave enough time to induce any potential 
nutrient stimulation. After 24 h, all bottles were spiked with ^^N2 
as described above, and incubation was continued under the same 
conditions for an additional 24 h. At the end of each incubation, 
the three treatments and control replicates were filtered as 
described above in order to measure N2 fixation rates, and 
amplification of the nifH gene (2010 only). Samples were also 
collected from bottles sacrificed at time zero in order to quantify 
background NOx and P04^~ concentrations at every station. NOx 
concentrations were also measured just after the NO3 additions 
in order to confirm the added concentrations at the beginning of 
the incubations (data not shown). 



PLOS ONE I www.plosone.org 



3 



December 2013 | Volume 8 | Issue 12 | e81265 



Aphotic N2 Fixation in the Pacific Ocean 



Mass spectrometric analyses 

The isotopic enrichment of particulate N after the incubation of 
seawater with ^^N2 was measured by continuous flow isotope ratio 
mass spectrometry of pelletized filters (Europa Integra-GN), 
calibrated every 1 0 samples using reference material (International 
Atomic Energy Agency [AIEA], Analytical Quality Control 
Services). The linearity of ^^N atom % as a function of increasing 
sample PON mass was verified as detailed in [40] on both 
natural and ^^N enriched material. This step is critical in ultra- 
oligotrophic environments or deep waters, where suspended PON 
concentrations are low. ^^N atom % was linear (Fisher test, 
p<0.01) between 0.20 and 39 |imol N, which is within the range 
of PON measured in all of our samples (0.27 to 4.91 |imole N 
depending on the station and depth). 

Detection and quantification limits for particulate N were 
calculated daily, as 3 times and 10 times the standard deviation of 
^^N analysis of blanks, respectively. Detection limits ranged from 
0.10 to 0.17 |imole N, and quantification limits ranged from 0.13 
to 0.26 |imole N, depending on the station. The ^^N isotope 
enrichment of a sample was calculated using the ^^N atom % 
excess over the ^^N atom % in samples taken from the same 
station at time zero, which was determined on bottles filtered 
immediately after adding ^^N2. We considered the results to be 
significant when ^^N excess enrichments were greater than 3 times 
the standard deviation obtained with ten AIEA references (^^N 
atom % >0.0005). The quantification limit of N2 fixation in this 
study was 0.01 nmol L ^ d \ If only one of the 3 replicate 
measurements was quantifiable, the average of the 3 replicates was 



forced equal to zero, in order to provide minimum estimates of N2 
fixation. 

In order to determine areal rates, N2 fixation measurements 
were trapezoidally depth-integrated from the summed products of 
the average of two adjacent rate measurements (including those 
equal to zero) with the depth interval between them. The standard 
deviation on the triplicates (2011 cruise) was also used for a 
trapezoidally depth-integration in order to obtain the standard 
deviation on integrated rates. 

Statistical analysis 

Controls and experimental nutrient treatments were compared 
using a 2-tailed non parametric Mann-Whitney mean comparison 
test (n = 3, 0( = 0.05, unpaired samples). 

Phylogenetic characterization of diazotrophs 

In order to characterize the potential diazotrophs present in the 
core of the OMZ that responded to the addition of glucose, nucleic 
acid samples were collected from triplicate bioassays during the 
2010 cruise for amplification of the nifH gene. At TO and at the 
termination of the experiment, bottles were immediately filtered as 
described in [41] onto 25-mm, 0.2-|im Supor filters (GE Osmotics, 
Minnetonka, MN), and immediately flash frozen in liquid N2. All 
filters were stored at — 80°C thereafter. 

DNA samples were extracted using the Qiagen All Prep kit 
(Valencia, CA), according to manufacturer's guidelines, with 
modifications to include freeze-thaw and bead-beating steps to 
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Figure 2. Horizontal and vertical distributions of hydrological and biogeochemical parameters during the 2010 cruise (R/V 
Atlantis). (A, B) temperature, (C, D) dissolved oxygen, (E, F) chlorophyll fluorescence in the northern transect (left panels) and southern transect 
(right panels). 
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disrupt the cells [42]. The wash steps of this protocol were 
automated using a QIAcube (Qiagen). DNA extracts were stored 
at — 20°C until use. 

Nested PGR amplification targeting a fragment of the nifH gene 
was carried out using degenerate primers nifHXA [43,44] using the 
reaction and thermocycling conditions described in [42]. Ampli- 
cons were purified using a QIAquick Gel Extraction Kit (Qiagen) 
and cloned using the TOPO TA Gloning Kit for Sequencing 
(Invitrogen, Garlsbad, GA) according to the manufacturer's 
guidelines. Purified recombinant plasmids containing partial nifH 
sequences were recovered from clones using the Montage Plasmid 
Miniprep96 Kit (Millipore, Billerica, MA) and sequenced using 
Sanger technology at the UG Berkeley DNA Sequencing Genter. 
All DNA extractions and as PGR preparations were performed in 
a PGR-amplicon free facility at UGSG described in [42]. 

Sequencher 5.1 sequence analysis software (Gene Godes 
Gorporation, Ann Arbor, MI) was used to remove vector 
contamination and low-quality reads from raw sequences. All 
resulting partial nifH sequences were imported into a curated nifH 
database (http://pmc.ucsc.edu/~wwwzehr/ research/ database/), 
translated into amino acid sequences, aligned to the existing 
hidden Markov model alignment using the Quick Align function, 
and nucleic acids were realigned to the aligned amino acids in the 
ARB software environment. Sequences generated from the 
nutrient addition bioassays were clustered at 97% nucleotide 
similarity using GD-HIT-EST [45]. Nucleic acid trees used the 
Jukes-Gandor correction for branch length. Trees generated in 
ARB were exported into iTOL for the display of associated 
metadata. All partial nifH sequences recovered were submitted to 
Genbank under Accession numbers KF5 15738 - KF5 15848. 



Results 

Hydrographic and nutrient profiles 

During both cruises, oceanographic conditions were consistent 
with active wind-driven upwelling off the coast of Northern Ghile 
and Peru (Fig. 1), associated with a vertically and horizontally 
extensive OMZ (Fig. 2G and 3B). During the 2010 cruise (R/V 
Atlantis, El Nino), the zone of decreasing dissolved oxygen 
(oxycline) was located at ca. 70-100 m at stations 11 and 9, and 
suboxic conditions ([O2] <20 |imol kg~^ [46]) were reached at 
125 m and 130 m respectively at these two stations. The suboxic 
zone expanded from 130-750 m at Station 11. During the 2011 
cruise (Fig. 3B, R/V Melville, La Nina conditions), the oxycline 
was shallower (ca. 30-40 m) at stations 1 1 and 9 compared to 
2010, and suboxic conditions were reached at 65 m at both 
stations. The suboxic zone expanded from 70-850 m at Station 
1 1. On the southern transect (Stations 1, 3 and 5, Fig. 2D and 4B), 
the water column was well oxygenated during both cruises (ca. 
200 |imol kg~ ^) over the first 200 m, and O2 concentrations 
decreased with depth to reach minimum values of ca. 
80 |imol kg~^ at 300 m depth. 

During the 2011 cruise, on the northern transect at Stations 7 
and 9 (Fig. 3D), surface NOx concentrations were ca. 6 |imol L \ 
increased quickly with depth (nitracline ca. 50 m) to reach ca. 
37 |imol in the core of the OMZ. Glose to the coast (Station 
1 1), the nitracline was shallower (ca. 25 m) than that of the oceanic 
stations, and NOx concentrations increased quickly to reach 
concentrations of ca. 30-40 |imol L ^ in the core of the OMZ. 
On the southern transect (Fig. 4D), NOx concentrations were ca. 
0.10 |imol L ^ over the first 100 m of the water column. The 
depth of the nitracline was 85 m and 130 m at Stations 1 and 5, 
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Figure 3. Horizontal and vertical distributions of hydrological and biogeochemical parameters during the 201 1 cruise (R/V Melville) 
- Northern transect (10°S). (A) temperature, (B) dissolved oxygen, (C) chlorophyll a fluorescence, (D) NOs" concentrations, (E) P04^~ 
concentrations, (F) Mean N2 fixation rates (n = 3). 
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Figure 4. Horizontal and vertical distributions of hydrological and biogeochemical parameters during the 201 1 cruise (R/V Melville) 
- Southern transect (20°S). (A) temperature, (B) dissolved oxygen, (C) chlorophyll a fluorescence, (D) NOa" concentrations, (E) P04^~ 
concentrations, (F) Mean N2 fixation rates (n = 3). 
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respectively, and concentrations increased progressively to reach 
values of ca. 40 |imol L~ ^ below the oxygen minimum and down 
to 2000 m. PO4 concentrations (Fig. 3E and 4E) followed the same 
trend as NOx in the northern transect, with surface concentrations 
of 0.40-0.60 |imol L ^ and a shallower phosphocline near the 
coast (Station 11) compared to open ocean stations. In the 
southern transect, surface PO4 concentrations were lower com- 
pared to those of the northern transect (0.03-0.30 |imol L ^) and 
the phosphocline was located deeper (ca. 150 m). 

During both cruises, chorophyll a fluorescence (Fig. 2E-F, 3C, 
4C) was highest at stations located along the northern transect. It 
was much lower in 2010 compared to 201 1, especially at stations 
nearest the coast of Peru on the Northern transect. 

Vertical profiles (0 to 2000 m) of N2 fixation 

During the 2010 cruise, N2 fixation was detected in 34 of 40 
samples representing all three stations (Fig. 5). The overall range of 
rates measured over the cruise was from the detection limit to 
0.80 nmol L ^ d \ The highest rates were measured in O2 
deficient waters at the oxyclines or in the core of the OMZ (Fig. 5), 
and reached values up to 0.57 nmol d~\ 0.6 nmol d~^ 
and 0.53 nmol d~^ at Stations 1, 11 and 9, respectively. 
Below the OMZ, rates were always measurable and were at 
1000 m depth 0.16 nmol L"^ d~\ 0.23 nmol L"^ d"^ and 
0.06 nmol L ^ d ^ at these stations. Integrated rates over the 
2000 m water column were 155 |imol N m ^ d ^ at Station 1, 
288 iimol N m"^ d"^ at Station 9, and 509 |imol N m"^ d"^ at 
Station 1 1 (Table 1). The average integrated rate over the cruise 
was 3 1 7 |imol N m ^ d \ Integrated N2 fixation rates in the 



aphotic zone accounted for 73 to 99% of the rates measured over 
the entire water column depending on the station. When 
considering all the stations, the average areal rate in the aphotic 
zone was 87% of the total rate over the entire water column 
(Table 1). 

During the 20 1 1 cruise, N2 fixation rates were significantly 
greater than zero in 140 of the 216 measurements made (Fig. 3F, 
Fig. 6). The overall range of rates measured was from detection 
limit to 0.26±0.12 nmol L ^ d \ In the northern transect 
(Fig. 3F, Fig. 6), the highest rates of N2 fixation over the vertical 
profiles were measured in the oxycline as in 2010, and mean rates 
(n = 3) reached 0.15±0.13 nmol L"^ d"^ at Station 7 and 
0.19±0.28 nmol L ^ d ^ at Station 9 at the oxycline. At station 
1 1 , the highest rates were found in surface waters 
(0.22±0.19 nmol d~^) but rates at the oxycline were also 
measurable (0.06±0.03 nmol L"^ d"^). Below the OMZ (ca. 
400-2000 m), rates were also measurable and ranged from 
0.00±0.01 to 0.21+0.13 nmol L"^ d~\ the highest rates being 
measured at station 9 at 1000 m depth. In the southern transect 
(Fig. 4F, Fig. 6), the rates ranged from 0.00 ±0.01 to 
0.26±0.12 nmol L ^ d \ the highest rates being observed at 
Station 1 just below the second oxycline at 750 m depth. At this 
station, aphotic rates were measurable at 450 m, 750 m, 1000 m 
and 1500 m depth. Integrated rates over the 2000 m water 
column ranged from 24±14|imolm ^ d ^ at Station 5 to 
118±87 iimol m"^ d"^ at Station 1 (Table 1). The average 
integrated rates over the 2011 cruise were 64 |imol m~^ d~\ 
Integrated N2 fixation rates over the aphotic zone accounted for 
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Figure 5. Vertical profiles (0-2000 m) of N2 fixation rates (nmol d~^) during the 2010 cruise (R/V Atlantis). Open circles: individual 
N2 fixation measurements at each depth. Black line: dissolved oxygen concentrations (|amol Kg~^) divided by 500 to fit on the same scale. 
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90% of total rates measured over the entire water column (Table 1) 
over the cruise. 

Nutrient addition bioassays in tlie core of tlie OMZ 

During the 2010 cruise, nutrient concentrations in the core of 
the OMZ (140 to 450 m) where experiments were performed 
ranged from 24.0 to 37.0 |imol L"^ for NO^ and 1.20 and 
3.00 |imol for P04^~ (data not shown). Mean N2 fixation 
rates in control bottles at Stations 5, 7 and 11 were 0.12±0.02, 
0.16±0.04 and 0.17±0.02 nmolNL"^ d"\ respectively (n = 3; 
Fig. 7). At the 2 most oceanic Stations 5 and 7, glucose 
amendments did not result in any significant increase of N2 



fixation (p>0.05). At Station 1 1 near the Peruvian coast, 
glucose amendments resulted in a significant (p<0.05) increase 
in N2 fixation rates by a factor of 3.2, to reach 0.56 ± 
0.04 nmol N L"^ d"^ (Fig. 7). 

During the 2011 cruise, nutrient concentrations in the core of 
the OMZ (320 to 475 m) where experiments were performed 
ranged between 33.12 and 38.72 |imol for NOx and from 
2.33 to 3.03 |imol L ^ for PO4. Mean N2 fixation rates in the 
control bottles ranged from 0.00 ±0.01 at Station 1 1 to 0.07 ±0.01 
and 0.07 ±0.04 nmol N L ^ d ^ at Stations 9 and 5, respectively 
(n = 3; Fig. 8). At Stations 1 and 9, N2 fixation rates were 
significantly (p<0.05) stimulated by simple carbohydrate additions 



Table 1. Area! N2 fixation rates (|imol N m ^ d ^) calculated from measurements performed in the euphotic and aphotic zones 
during the 2010 and 2011 cruises. 



Station 


Total 0-2000 m (|imol N d^) 


Euphotic (jimol N m ^ d ^) 


Aphotic (jimoi N m ^ d ^) 


2010 


Station 9 


288 


49 


239 


Station 11 


509 


8 


501 


Station 1 


155 


71 


84 


Average 


317 


43 (13%) 


275 (87%) 


2011 


Station 7 


50±39 


10±10 


40 ±29 


Station 9 


101±79 


09 ±04 


92 ±74 


Station 11 


25±17 


04 ±02 


21±14 


Station 5 


24±14 


05 ±05 


19±09 


Station 1 


118±87 


01 ±01 


117±86 


Average 


64 


06 (10%) 


58 (90%) 



Uncertainties are derived from tlie standard errors of triplicate measurements for the 2011 cruise. 
doi:1 0.1 371/journal.pone.0081 265.t001 
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by a factor of 5.5 and 4.6, to reach 0.14±0.07 and 
0.30±0.30 nmol N L"^ d~\ respectively (Fig. 8). At all stations, 
the addition of amino acids resulted in a significant (p<0.05) 
increase in N2 fixation rates, by a factor of 4 to 7. The highest 
rates were reached at Station 9 after AA additions with 
0.27±0.08 nmol N L"^ d"^ ATP addition never resulted in 
any significant increase of N2 fixation rates (p>0.05) and NOs" 
additions never resulted in any decrease of N2 fixation rates 
(p>0.05). However, at Station 7, NOs" additions resulted in a 
significant (p<0.05) increase of N2 fixation by a factor of 7 to 
0.07±0.02 nmol N L"^ d~\ 

Phylogenetic characterization of diazotrophs in 2010 
glucose addition bioassays 

A full phylogenetic characterization of diazotrophs in the upper 
200 m of the ETSP water column was performed during the same 
cruises and is detailed in a companion paper [42] . In this study we 
report the complementary phylogenetic characterization of 
samples from the core of the OMZ (Fig. 8). Partial ni/H sequences 
recovered during deep glucose addition bioassays during 2010 at 
Stations 5, 7 and 1 1, indicated that diazotrophs were present in the 



deep waters of the OMZ. The diazotrophic community was 
comprised of non-cyanobacterial diazotrophs affiliated with nifH 
clusters IK (predominantly comprised of cx-proteobacteria), IG 
(predominantly comprised of y-proteobacteria), and 3 (sulfate 
reducing genera of the 5-proteobacteria as well Clostridium spp.. 
Vibrio spp, etc.) (Fig. 9). Clear differences exist between OMZ 
diazotrophic community composition at each station. The Station 
5 community was dominated by nifH cluster IK sequences, many 
of which are closely related to a phylotype (94-97% nucleic 
acid similarity) originally reported at Hydrostation S (North 
Atlantic) from a depth of 1000 m (BT5167A10 (DQ481253) [47]), 
although a few putative y-proteobacteral (IG) sequences were also 
recovered that affiliated with yETSPS, a cluster recovered from 
the ETSP [42]. Although the lowest number of total sequences was 
recovered from Station 7, they were mainly affiliated with cluster 
IG, along with a few IK sequences. In contrast, clone libraries 
from Station 1 1 were dominated by cluster 3 sequences, along with 
a few IG sequences, but no IK sequences (Fig. 9). 

Despite being prevalent in clone libraries, both (xETSPl and 
cIITETSP groups were not detected at abundances great enough 
to quantify using Taqman® qPCR assays during the course of 
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Figure 7. Effect of glucose additions on mean N2 fixation rates 
(n = 3) during the Atlantis cruise (2010). Data from bioassay 
experiments performed in the OMZ at Stations 5, 7 and 11. The error 
bars represent the standard deviation of triplicate incubations. 
Treatment means were compared using the 2-tailed non parametric 
Mann-Whitney mean comparison test (n = 3, a = 0.05, unpaired sam- 
ples). Means that are significantly different (p<0.05) from the control 
are labeled with an asterisk. 
doi:1 0.1 371/journal.pone.0081 265.g007 



these experiments (see Figuer SI in File SI). Because the 
abundances of these targets did not increase as a result of nutrient 
amendments, it is difficult to speculate whether any of them were 
responsible for the increased N2 fixation rates we measured after 
glucose addition. 

Discussion 

Active N2 fixation in deep and NOs-rich waters of the 
ETSP 

In this study, we measured during 2 consecutive years N2 
fixation in surface waters affected by the OMZ, but reveal that N2 
fixation below the euphotic zone is more important: 87 and 90% 
of total areal N2 fixation were measured in the aphotic zone in 

2010 and 201 1, respectively. In these aphotic layers, volumetric N2 
fixation rates were relatively low (<1.00 nmol N d~^), but 
when integrated over the whole aphotic layer, they ranged from 
84 to 501 iimol N m"2 d"^ in 2010 and from 19±09 to 
117±86 iLimol N m"2 d"^ in 2011 (Table 1). In 2011, rate 
measurements were replicated (triplicates) and calculations 
performed very carefully using a real TO for every depth. These 

201 1 measurements are thus more reliable than those measured in 
2010. These measurements in aphotic waters add new information 
compared previously published studies [26,28] in the area. The 
hypotheses explaining the persistence of N2 fixation in these high 
NOx (ca. 40 |imol L ^) environments are largely developed in the 
companion paper [28]. First, fixed N loss processes occur in this 
region [31,48], creating a deficit of N relative to P, which is 
potentially favorable for N2 fixation [24] . In particular, anammox 
removes NH4^, which has an immediate inhibitory effect on N2 
fixation [49]. Secondly, N2 fixation is an anaerobic process [50] 
due to the irreversible inactivation of the nitrogenase enzyme by 
O2 [51]. It is possible that the low O2 concentrations in the OMZ 
and down to 2000 m contribute to the protection of nitrogenase 
[52], decrease the energy cost of maintaining intracellular 
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Figure 8. Effects of nutrient additions on mean N2 fixation rates (n = 3) during the Melville cruise (2011) (C: Control, AA: Amino 
acids, CH: Carbohydrates). Data from bioassay experiments performed in the OMZ at (A) Stations 7, (B) Station 9, (C) Station 1 1, (D) Station 5, (E) 
Station 1, and (F) Station 6 during the Melville cruise (2011). The error bars represent the standard deviation of triplicate incubations. Treatment 
means were compared using the 2-tailed non parametric Mann-Whitney mean comparison test (n = 3, a = 0.05, unpaired samples). Means that were 
significantly different (p<0.05) from the control are labeled with an asterisk. 
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Figure 9. Results from the molecular analyses of the 2010 glucose amendment experiments. Neiglibor joining tree of partial nifH 
nucleotide sequences. Nodes are labeled with nifH cluster designations according to the convention established [74]. The number of sequences 
recovered fronn stations and treatnnents for each phylotype are indicated in boxes to the right of the tree. 
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anaerobiosis [53], and thus facilitate N2 fixation. Finally, redox 
conditions in the OMZ favor the equilibrium formation of the 
most bioavailable form of iron Fe^^ [54], which could help to 
support the high Fe requirements of nitrogenase [16,17]. For these 
reasons, OMZs and deeper waters may represent favorable 
ecological niches for N2 fixation, as shown in this study. 

Potential impacts on N budgets in the ETSP 

Aphotic N2 fixation is currently ignored in oceanic N budgets 
based on biogeochemical rate measurements. However, this 
dataset indicates that rates in aphotic waters of the ETSP are of 
the same order of magnitude than those commonly measured in 
the tropical and sub-tropical NO 3 -depleted surface ocean 
(Table 2), where N2 fixation has commonly been studied. The 
potential significance of the N2 fixation rates measured in our 
study can be evaluated by comparing them with fixed N losses via 
denitrification and anammox measured in the same region. N 
losses in the ETSP have been estimated to range from 9 to 25 Tg 
N yr"' (Table 3, [31,55,56,57]) in the upwelling area extending 
175 km offshore, and 1860 km along the Peruvian-Chilean coast 
with an area extant of 3.26x10^^ m^ [31,56]. If we consider the 
same spatial extent for N2 fixation, this process could potentially 
add 0.04 to 0.9 Tg N yr"^ (Table 3), counterbalancing 0.16 to 
10% of the estimated N loss processes in this area (calculations 
have been performed only using numbers from the 2011). 
However, the anammox and denitrification measurements men- 
tioned above [31] were performed under conditions of excess 



substrate availability and therefore represent maximum estimates 
of N loss rates. In contrast, the ^^N2 bubble method used to 
quantify N2 fixation [58] may underestimate rates [38,39]. 

Secondly, denitrification and anammox are restricted to 
subsurface suboxic or anoxic waters [59], whereas N2 fixation is 
not. Further, denitrification and anammox appear to be restricted 
to the coastal upwelling system within ca. 175 km of the Peruvian- 
Chilean coast (the few data available at open ocean stations 
indicate that N loss processes were below detection limit during the 
2010 cruise, Hamersley et al., (Pers. Com.)). N2 fixation in the 
ETSP is active over a much greater spatial extent than N 
loss processes. If we consider the spatial extent of the N2 
fixation measurements in the ETSP covered by our cruises 
(2.23 X 10^^ m^), we estimate that N2 fixation could potentially add 
0.3 to 1 Tg N to the system in this area and therefore could 
compensate for up to 11 % of the estimated N loss processes in the 
upwelling region of the ETSP (Table 3) (without taking into 
account methodological under- or overestimations). These esti- 
mates of N gains are the minimum ones calculated by taking into 
account only the 2011 cruise. If we take into account the 2010 
cruise, N2 fixation could potentially compensate up to 0.3 to 7 Tg 
N to the system (Table 3) (i.e. up to 78% of N losses). N2 fixation in 
deep waters of the ETSP may be a significant source of N for the 
ETSP, and needs to be taken into account in future N budgets. 
Further coupled measurements between N gain and loss processes 
at the same stations/depths need to be performed to better 
constrain the magnitude of N gains in this region. 
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Table 2. Examples of published studies showing the range of oceanic N2 fixation areal rates measured in some contrasting 
oceanic environments. 



Location 


Areal rates (jimol m ^ d ^) 


Integration depth (m) 


Reference 


Hypoxic basin (Southern California Bight) 


150 


885 


[37] 


ETSP coastal OMZ 


7-190 


120 


[26] 


ETSP 


0-148 


150-200 


[28] 


ETSP subtropical gyre 


12-190 


150-200 


[14] 


Eastern North Pacific gyre 


520 


mixed layer 


[20] 


North Atlantic 


59-898 


15 {Trichodesmium bloom) 


[3] 


ETSP aphotic zone 


19-501 


2000 


This study 



doi:1 0.1 371/journal.pone.0081 265.t002 



Effects of nutrients on N2 fixation 

The diazotrophic community of the ETSP characterized in this 
study, as well as in a companion study [42] is comprised of an 
assemblage of non-cyanobacterial diazotrophs, and little can be 
inferred about their metabolism from partial nifH sequences. 
However, we performed nutrient addition bioassays using 
molecules representing common labile components of the 
dissolved organic matter pool in marine waters (simple carbohy- 
drates, amino-acids and ATP), which shed some light on nutrient 
control of N2 fixation in the core of the OMZ. Our results 
indicated that simple carbohydrate additions significantly stimu- 
lated N2 fixation at stations located nearest the coast during both 
cruises and at Station 9 during the 2011 cruise (Figs. 7, 8). In 
OMZs, organic C is largely supplied by vertical flux of planktonic 
production from shallower layers or by horizontal transport [60] . 
Thus this supply of organic G is not constant but rather episodic, 
which could explain why N2 fixation appears so variable in space, 
in time, and between cruises and years, as reported in the present 
study and by [26]. This seems to be the case for N loss processes as 
well, since organic C supply has been correlated with regional 
denitrification [60] and anammox [31,61] rates in OMZs. Ward et 
al. [62] demonstrated that denitrification rates were significantly 
stimulated in the OMZ of the ETSP by organic C additions. To 
our knowledge our study is the first designed to study the response 
of diazotrophs to nutrient additions in the OMZ. In surface 
waters, significant stimulation of N2 fixation rates by glucose 
additions have been reported during the same cruise at Station 9 
[28]. A significant stimulation of bacterial production after glucose 
amendments in surface waters of the Chilean up welling system 
have also been reported [63]. Finally, in surface waters of the 
southwest Pacific [64] , reported a significant increase of nifH gene 
copies of unicellular diazotrophic cyanobacteria such as Group A 
(UGYN-A) and Crocosphaera after glucose and mannitol additions, 
hypothesizing that this capacity may allow conservation of energy 



by rapid uptake and recycling of sugars. However, it has to be 
noted that the large variability in the response to carbohydrates 
addition (high standard deviation at Station 1 1 for example) could 
be explained by the fact that it may be coincidental whether the 
taxa that benefit from the enrichment possess the nifH gene. 

Because the organic C molecules tested here are also energy- 
rich molecules easily entering catabolic pathways, one could 
interpret our results to be indicative of limitation either by energy 
or by assimilative C availability. However, in our experiments N2 
fixation was not stimulated by ATP additions at any station, 
indicating that C and not energy might have been the proximate 
limiting factor. In some oligotrophic P-limited environments, ATP 
is also a source of P for bacteria and uptake rates of ATP exceed 
those of glucose [65]; however, in OMZs, P is not limiting relative 
to N, which may further restrict the ability of ATP to stimulate N2 
fixation rates in our bioassays. In contrast, the addition of free 
amino acids stimulated N2 fixation at all stations tested; this has 
also been shown in aphotic oxynenated waters of the Red Sea 
[66] . Amino acids are a source of both C and N, and it has been 
suggested that it is energetically advantageous for microbes to use 
preformed compounds such as amino acids rather than glucose as 
C sources [67]. In terrestrial legume-rhizobium symbioses, the 
diazotrophic bacteria assimilate amino acids such as glutamic acid 
provided by the host, which facilitate both dicarboxylate oxidation 
and ammonium assimilation into asparagine [68]. In Azospirillum 
sp.^ additions of glutamic acid also stimulated N2 fixation activity 
[69] by providing a C and energy source to the diazotrophs, while 
N was still provided via N2 fixation. It may be that similar nutrient 
assimilation dynamics are occurring in diazotrophs in the ETSP 
OMZ. 

Ambient NOx concentrations were high (ca. 30-40 |imol L ^) 
at all stations where nutrient additions were performed, and NO3 
additions (8 |imol L ^) never resulted in N2 fixation inhibition at 
any station. As the metabolic potential of diazotrophs present in 



Table 3. Comparison between estimated fixed N losses via denitrification and anammox and fixed N gains via N2 fixation 
(estimated from the 0-2000 m depth integrated rates measured in this study) in the ETSP. 



Area considered N losses N gains (based on 2011 cruise) N gains (based on 2010+2011 cruise) 



(m^) 


(Tg.yr^) 


(Tg.yr^) 


(Tg.yr^) 


^3.26x10" 


^^9-25 


0.04-0.9 


0.05-1.1 


2.23x10^^ 


9-25 


0.3-1 


0.3-7 



^Upwelling area extending 175 km from, and 1860 km along the Peruvian-Chilean coast [31,56]. 

^^Estimates from [31,55,56,57]. 

doi:1 0.1 371/journal.pone.0081 265.t003 
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the OMZ have not yet been fully characterized, we do not know if 
they possess genes for reduction and assimilation of NO 3 or 
NO2 . Detailed studies at the single cell level would be needed to 
characterize the metabolism of these organisms and understand 
why microbes fix N2 in the presence of so much NO 3 . In 
addition to possible energy and C, N and P sources derived from 
molecules like amino acids, carbohydrates or ATP, electron 
sources and donors are also very important to know for 
characterizing the physiology of the diazotrophs present in the 
ETSP. Molecules like O2, NOs" and less favorably S04^~ are 
common electron acceptor and they are used for different types of 
respirations like aerobic respiration, or anaerobic denitrification 
and sulfate reduction. These different respiratory pathways 
potentially supporting N2 fixation are performed by organisms 
with different physiology which each have their own environmen- 
tal sensitivities for ffxing N2. 

Phylogenetic characterization of diazotrophs 

Our characterization of the diazotrophic community in the core 
of the OMZ revealed the presence of potential N2-fixing 
heterotrophs based on the presence of the nifH gene. We did not 
detect the cyanobacterial diazotrophs commonly found in other 
regions of the open ocean; in contrast, most of the niJH genes 
amplified from the OMZ clustered with a-, y- and 5-proteobac- 
teria. This result is consistent with the observations of Turk-Kubo 
[42] in the upper 200 m of the ETSP water column, where 96% of 
sequences were also affiliated with proteobacteria. Based on these 
results, and other studies conducted in the ETSP and the South 
Pacific Gyre [14,26,70], it is clear that the ETSP diazotrophic 
community is different from other well-studied tropical and sub- 
tropical oceans such as that of the North Pacific, North Atlantic 
and Indian Oceans. The cyanobacterial diazotrophic phylotypes 
commonly found at high abundances in these other ocean 
provinces appear to be either sporadically present at low 
abundances (i.e. Trichodesmium, UCYN-A), or undetected altogeth- 
er (i.e. UCYN-B, diatom-diazotroph associations) in the ETSP. 

The amplification of diverse non-cyanobacterial ^z/?/-containing 
organisms from OMZ waters in the ETSP affiliated with nifH 
clusters IK, IG and 3, is consistent with the findings of other 
studies conducted in anaerobic waters [26,36,71] and in 
abyssopelagic waters [47]. However, the results from this study 
underscore the difficulty inherent in identifying the diazotrophic 
community responsible for N2 fixation rates. It is important to note 
that nifH cluster IK sequences have been reported as contam- 
inants in many studies, including a study in the ETSP [42]. 
However, none of the sequences recovered here had greater than 
90% amino acid similarity and 83% nucleic acid similarity to 
reported contaminants. Nevertheless, as a result of the use of 
highly degenerate primers and nested PGR cycles necessary to 
amplify this important but low-abundance gene target, contam- 
ination must always be considered as a source for heterotrophic 
diazotroph sequences, whether from PGR and DNA extraction 
reagents or from sampling or handling procedures, despite the 
screening of PGR and reagent blank controls as in this study. 

Furthermore, although it is clear that a diverse assemblage of 
non-cyanobacterial /2z/?/-containing organisms are present in the 
OMZ of the ETSP, the best methodologies currently available to 
characterize dominant members of the diazotrophic community 
(PGR amplification using degenerate nifH primers) often identify 
organisms present at extremely low levels when targeted using 
quantitative approaches (i.e. qPGR; Fig. SI in File SI) [42,47,71]. 
This, in turn, makes it difficult to argue that these organisms are 
capable of fixing N2 at cell-specific rates great enough to account 
for measured bulk rates. An analysis of the expected N2 fixation 



rates based on abundances and plausible cell-specific N2 fixation 
rates in the ETSP discussed in [42] indicate that these 
proteobacteria are unlikely to be responsible for all the measured 
bulk rates and therefore other N2-fixing organisms could be 
responsible for a part of N2 fixation in this region but may remain 
uncharacterized. Identifying which organisms are actively tran- 
scribing nifH using techniques such as reverse transcription (RT)- 
qPGR might provide more insight into which diazotrophic taxa 
are actively fixing nitrogen. However, the challenge of identifying 
which organisms are important N2-fixers remains the same when 
designing qPGR primers from sequences derived from RT-PGR 
based clone libraries, and are further convoluted by potentially low 
transcript abundances per cell and/ or the timing of sampling with 
respect to diel changes in nifH expression (even in the case of 
heterotrophs). 

Conclusions 

This study provides one of the first estimates of N2 fixation rates 
in aphotic waters of the ETSP. It reveals that N2 fixation in 
aphotic environments is the largest contributor to total areal N2 
fixation in ETSP. N2 fixation in high [NO 3 ] environments 
remains an enigma as it requires an additional energetic cost 
relative to NO 3 or NH4 . Further physiological studies are 
needed to understand the physiological regulation of N2 fixation, 
especially on newly discovered diazotrophic organisms. Gontrary 
to N2 fixation performed in euphotic layer which sustains new 
primary production [3], aphotic N2 fixation may sustain organic 
matter remineralization. These new sources of N could potentially 
compensate for as much as 78% of the estimated N loss processes 
in ETSP, indicating that they need to be taken into account in 
marine N budgets. Phylogenetic studies confirm the presence of 
diazotrophs in the deep waters on the OMZ, which are distinct 
from cyanobacterial phylotypes commonly found in surface 
oligo trophic waters of the tropical ocean. Organic and inorganic 
nutrient addition bioassays reveal that amino acids and simple 
carbohydrates stimulate N2 fixation in the core of the OMZ, and 
the episodic supply of these nutrients from upper layers may 
explain the large temporal and spatial variability of N2 fixation in 
the ETSP. Research on marine heterotrophic N2 fixation is at its 
beginning and significant progress needs to be made in the 
refinement of the methods to estimate planktonic N2 fixation in 
OMZs (^^N2 bubble method versus N2-enriched seawater) from 
bulk measurements to single cells analysis. The ^^N-enriched 
seawater method should be coupled to oxygen-free and trace 
metal-clean procedures to provide more accurate estimates. 
Progress also needs to be made in the characterization of the 
community responsible for N2 fixation in these deep waters, as well 
as the control of their population dynamics by the supply of 
organic matter. Estimates of global N2 fixation based on field 
measurements [5,72] are presently lower than geochemically- 
based (nutrient stoichiometry and isotopic ratio) estimates [73]. 
Taking into account deep N2 fixation might help to resolve some 
of this discrepancy. However, progress also needs to be made in 
the quantification of N loss processes, as recent studies indicate 
that they may be less sensitive to oxygen than previously thought 
[59], further complicating the N budget in the ETSP. In future 
studies, N gain and loss measurements need to be coupled in space 
and time to further resolve the N budget in the ETSP. 
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